Quasiperiodic dynamics of coherent diffusion: a quantum walk approach 
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We study the dynamics of a generalization of quantum coin walk on the line which is a natural 
model for a diffusion modified by quantum or interference effects. In particular, our results provide 
surprisingly simple explanations to phenomena observed by Bouwmeester et al. (Phys. Rev. A 61, 
13410 (1999)) in their optical Galton board experiment, and a description of a stroboscopic quantum 
walks given by Buershaper and Burnett ( quant- ph/0406039 1 through numerical simulations. We 
also provide heuristic explanations for the behavior of our model which show, in particular, that its 
dynamics can be viewed as a discrete version of Bloch oscillations. 
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The behavior of quantum analogs of classical systems 
with chaotic dynamics has recently attracted much at- 
tention of both theoretical and experimental physicists. 
One of the most widely studied non-classical phenomena 
emerging in such systems are quantum suppression (dy- 
namical localization) and enhancement (quantum reso- 
nance) predicted for various models of spectral diffusion. 
Both suppression and enhancement have been observed 
for a model of <5-kicked rotor in an experiment in which 
cold atoms interact with a pulsed optical lattice [1-4] 
and, recently, in superconducting nanocircuits Q. An- 
other approach to quantum chaos have been presented by 
Bouwmeester et al. |6|]. They realized an optical version 
of a Galton board (OGB), and studied how the diffusion 
in OGB is influenced by a choice of two parameters - the 
diabatic transition probability D and the relative phase 
$ (precise definitions of D and <E> are given below) . The 
dynamics of OGB turned out to be very different from a 
classical Galton board which corresponds to a standard 
random walk and results in Gaussian diffusion. In partic- 
ular, the authors of the OGB experiment concluded that 
"suppression of diffusion can occur in the form of almost 
perfect recurrences of the initial level population" 
Such an effect had been already predicted by Harmin 
in the adiabatic limit (D <§; 1), but, since these recur- 
rences must vanish in the strictly diabatic limit (D = 1), 
it came as a surprise that they were observed even for 
the transition probability D as high as 0.8. Prompted by 
these observations the authors speculated that the sys- 
tem has two different modes of behavior, adiabatic and 
diabatic. A similar picture of the dynamic of OGB was 
presented by Torma Q who even provided some heuris- 
tic arguments suggesting that, in fact, the transition be- 
tween adiabatic and diabatic phases may be similar to the 
phase transition in two-dimensional Ising model. We also 
remark that in the OGB experiment the diffusion is sup- 
pressed for rational values of the relative phase 0/27T, for 
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which one would rather expect an enhancement caused 
by quantum resonance effects. 

In this paper we model coherent diffusion using a gen- 
eralized coin quantum walks (GCQW) providing a consis- 
tent, although somewhat unexpected, explanation for re- 
currences phenomena observed in the OGB experiment. 
In particular, our model predicts that the recurrences 
changes with the transition probability D in a smooth, 
continuous way which refutes the conjecture that in the 
dynamics of such a system one can distinguish two dif- 
ferent adiabatic and diabatic modes of behavior. More- 
over, our approach explains nicely why and when we can 
observe recurrences and ballistic diffusion; it also pre- 
dicts some new multiple recurrences effects which, most 
likely, can be verified experimentally within current tech- 
nology. Besides of the quantitative analysis we offer a 
simple intuitive explanation of basic features of coherent 
diffusion dynamics which explores similarities between 
GCQW and Bloch dynamics. Let us also remark that 
the quantum walk framework and results presented here 
are not restricted to particular implementation of coher- 
ent diffusion such as OGB. In fact, we believe that our 
approach can capture many general properties of diffu- 
sion modified by quantum or interference effects. 

The optical Galton board (OGB) studied by 
Bouwmeester et al. [|| consisted of a ladder of equally 
spaced levels periodically coupled to each other via 
Landau-Zener crossings. Due to the structure of equally 
spaced levels OGB resembles a <5-kicked harmonic oscil- 
lator. However, as was pointed by Bouwmeester et al. 
[6(, it can be assumed that Landau-Zener crossings in- 
duce transitions between neighboring levels only. Thus, 
the OGB experiment can be as well considered as an im- 
plementation of a coined quantum walk 9] . Let us recall 
that a coined quantum walk (CQW) is a quantum analog 
of a random walk on graph proposed by Aharonov, Davi- 
dovich and Zagury |10|, extensively studied for the last 
few years in the hope of constructing efficient quantum 
algorithms (see 01 )• Buerschaper and Burnett have 
noticed that CQW can also be useful for modeling and 
studying quantum chaos. Here we follow this approach 
and characterize the coherent diffusion in the terms of 
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generalized CQW. In a model of CQW on the line 01 
the evolution of the system is determined by a repeat- 
ing action of a unitary operator U defined on a tensor 
product of two Hilbert spaces H c ® H n . The base of H c 
consists of the coin states |c), (c = 0, 1), whereas H n is 
spanned by vectors \n), (n € Z), which correspond to 
the position of a particle on the line. In a single step of 
the evolution of the system we toss a coin, i.e., change 
the coin state |c), and move the particle to one of the 
two neighboring state of the line determined by |c). Let 
u = (E c l c >( c l ® S C )(C I), where C € 5C/(2) is an 
arbitrary coin tossing operator, and S c is a translation 
operator defined as S c \n) = \n + (— l) c ). The state vec- 
tor |\&) of the system evolves in (discrete) time t accord- 
ing to |*(t)) = f/*|*(0)>. Note that if the state of the 
coin is measured after each step, our model corresponds 
precisely to the classical random walk. 

If a quantum walk is to describe OGB it must take 
into account the phase relations, so the above mentioned 
basic model has to be generalized. This can be done by 
replacing U by = (£ c \c)(c\ <g> S ccj> )(C ® I), where 
•Sc^l"-) = e l ^^\n + (— l) c ), and 4>(n) is a phase acquired 
by the particle at the position n. We remark that the 
generalized coin quantum walk (GCQW) defined above 
is consistent with a recently proposed stroboscopic quan- 
tum walk (SQW) 12]. In the SQW model a particle 
walks on a space consisting of eigenstates of some Hamil- 
tonian H, and the evolution is periodically perturbed 
with an operator U. Both SQW and CQW models can 
be obtained from GCQW by taking <fr(n) = (n\H\n)t p 
(where t p is a period of perturbation) , and <fi(n) = const, 
respectively. Below we concentrate on harmonic case 
4>(n) — n$, where, furthermore, $ is rational fraction 
of 2ir, i.e., $ = 2ir^, with coprime q and p (although we 
shall make some assertions about the irrational case as 
well) . In order to study GCQW we assume that a parti- 
cle walks not on the line but on a long cycle of length N, 
where N is a large multiplicity of p so that 4>(N) = </>(0). 
Thus, here and below n € Zjv and addition in the def- 
inition of Serf, is taken modulo N . Obviously, for times 
t < N/2, the behavior of both the line and the iV-cycle 
models is identical. 

Let a coin tossing operator be in the form C = 
Y, c ((-l) c rf|c) (c| + a\c) (1 - c|) (if D is the diabatic tran- 
sition probability in the OGB experiment, then d = \f~D 
and a 2 + d 2 = 1). For the eigenvalues rj us , (j — 
0,1,..., N/p- 1; u = 0,1; s = 0,1,..., p - 1), of the 
unitary operator U$, for the ./V-cycle model we have got 
r jus = L0pZj u for an odd p, and r jus = ujp(2iXjZ ju - 1) for 

an even p, where tu p = e 27Tl / p , Zj U = ((— l) u ^Jl — X 2 — 

iXj^ 1 ^ , Xj — (iofrsin ^S^, the effective amplitude is given 
by d e f[ = d T / 2 , and T — p for an even p while T = 2p if 
p is odd. 

The probability that the particle returns to its ini- 
tial state is given by P(t) = |(*(t)|*(0))| 2 . Let 
f(t) = (*(i)|*(Q)) = Ej^^usv^, where A jus = 
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FIG. 1: The recurrence probability P(T) as a function of 
the diabatic transition probability D for p = 9. The curve is 
given by (|5J, the dots are found by numerical simulations. 

| (jus\ty(0)) \ 2 , and \jus) is the eigenvector corresponding 
to the eigenvalue rj US . 

In the adiabatic limit, when d — > 0, we have rj us = ±1 
independently of j, u, and s, so in this case P(T) = 1 as 
predicted by Harmin |2|. For d > we have 

f(T) = A+Re(Qu) + ]T Ajlm(Q u ) (1) 

j j 

where ( ju = rj us = ±(l - 2X 2 - 2^(-l) M A J yjl - X 2 ) (the 

sign depends on parity of p), and — s (Ajo s ± Aji s ) . 
Due to symmetries of eigenvectors second sum in Q van- 
ishes, whereas A^ = p/N for every j = 0,1, . . . , N/p — 1. 
Consequently, 

P(T) = (l-£ cff ) 2 , (2) 

where the effective probability of the diabatic transition 
is defined as 

D cS d 2 cff = D T ' 2 . (3) 

Note that, contrary to the suggestions that in the dy- 
namic of OGB there exist two sharply distinguished adi- 
abatic and diabatic phases 0,0, P(T) is a smooth func- 
tion of D (see Fig. [TJ. In fact, for every given D, one 
can always take p large enough so that P(T) = 1 and 
obtain an almost perfect recurrence, although, perhaps, 
due to decoherence effects, for such a p this phenomenon 
can hardly be verified experimentally. 

In order to study the behavior of GCQW for large t, 
let us calculate the probabilities P(kT) of multiple recur- 
rences. Approximating C^ u (up to the second order term) 

by Cj U = ±exp(2j(-l) u fcAj) leads to 

P(kT) = J 2 (2k^/D^), (4) 

where Jq{x) is a Bessel function (Fig. . Thus, the prob- 
ability P(kT) does not decrease monotonically with D e g 
but oscillates with this parameter, converging to as 
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FIG. 2: The probability of a multiple recurrence P(kT) as a 
function of k for D = 1/2 and p = 10. The curve is given by 
@l , the dots are found by numerical simulations. 

k — > oo. Since the first zero of Jq(x) occurs at x ~ 2.4, 
the decay of recurrences can be roughly characterized by 
time 

r(p)~1.2T£>- 1/2 . (5) 

This fact explains nicely numerical simulations for D = 
1/2 given by Buerschaper and Burnett 12] who, instead 
of recurrences, analysed the behavior of the related char- 
acteristic of a random walk: its standard deviation a. For 
example r(10) ~ 67.8, while r(20) - 768. Hence, the os- 
cillations of a for p = 10 die quickly (Fig. 2b in [l^l. 
while for p = 20 they can still be seen for t ~ 200 (see 
Fig. 2c in 01). In Fig.|31we show the striking difference 
in a dynamics for odd and even values of p. In the even 
case, when D e f{ = D p / 2 , (e.g., r(16) ~ 307.2), oscillations 
of a are damped more efficiently then in the odd case in 
which D e g = D p (t(15) ~ 6516.7). Numerical simula- 
tions suggest that for t ^S> r(p) a ballistic diffusion occurs 
with a proportional to t. In fact, the system behaves 
similarly to the basic model of CQW with D replaced 
by D c g; for instance, the formula a = ty 1 — \/l — D of 
Konno [1J| in the case $ = 2ir^ becomes 

cr = ty / l- v r l-D cS . (6) 

This formula holds for both odd and even p, although, as 
we have just observed, because of difference in D e ff, the 
value of a/t for p and p+ 1 can be dramatically different. 

One way to explain the behavior of GCQW described 
above is to compare this system to one with a ID peri- 
odic potential. The nodes of the TV-cycle correspond to 
potential wells. Since we use two sided coin in the corre- 
sponding system we have two bands of Bloch levels, each 
consisting of N levels, and the width of each band pro- 
portional to the tunnelling amplitude d. For $ = the 
symmetry of the system is given by Z^r. If $ = 27r^ the 
symmetry group reduces to Z/vy p and one can view such 
a system as a family of N/p clusters of wells, each with 2p 
energy levels, and the tunnelling amplitude d e g = d p , (so, 
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FIG. 3: The evolution of the standard deviation a for q = 1, 
p = 16 (thin line) and p = 15 (thick line). 
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FIG. 4: The levels of quasi-energy as a function of D e ff for 
p = 3, iV = 21. 



the width of the new band is proportional to d c g) . The 
quasi-energies Ej US of the system are given by Arg(rj US ) 
and, since the system evolves in discrete time, the whole 
spectrum of quasi-energies is contained within range of 
27r. Thus, at deft = 1, the width of a single band is 
2n/2p. Consequently, for d c g < 1, the width of each band 
should be approximately ^d c s- Our analytical compu- 
tations agree with this heuristic perfectly (Fig.^J. This 
picture is valid only for p odd; when p is even an addi- 
tional degeneration emerges and the resulting system has 
only p quasi-energy bands; so, in this case, d e g = 
In both cases, if d c s is small, the energy bands are very 
narrow, the resulting system is nearly harmonic and a 
periodic behaviour is expected. On the other hand, due 
to the uncertainty relation, the non-harmonic effects are 
expected to emerge only after time t larger than the reci- 
procity of the band width. Note that it is consistent with 
r(p) given by (J5J. 

Let us remark that although we study the case $ = 
none of the results presented so far have depended 
on q. The reason for that is that we have considered 
only " almost perfect" recurrences, considering only times 
t which are multiplicities of T. However, one can also 
observe "imperfect" recurrences for which the probabil- 
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ity of recurrence P(t) has a visible maximum, yet typ- 
ically is smaller than for the perfect ones. The exis- 
tence of such an imperfect occurrences is caused by the 
fact that our system is discrete, whereas in the con- 
tinuous version of the problem the maximum of P(t) 
might occur at t which is not an integer. From our 
previous considerations it comes as no surprise that an 
excellent candidate for a continuous version of our dy- 
namics is a Bloch system. To make this analogy more 
transparent, consider an equivalent, more symmetric, 
coin operator C — ^ c (id\c){c\ + a\c)(l — c|), and let 
= J2c n a c,n{t)\c)\n) . Then the evolution equation 
\^ (t)) — (0)} leads to the recursive equation 

Oc,„(t + 1) - a c , n (t - 1) = (e^ 2 "- 1 )* - l)a c , n (t - 1) 
+ ? :de^ 1 '*(a c ,, i „ 1 (i) + a C: „ +1 (<)) , 

for c = 0,1. Note that in the case of $ = an analogous 
equation has been studied in [HHE]- If n <C f = 2n/<& 
this equation can be approximated by the differential 
equation 



da c 



df 



in<fra c 



id 



(oe c 



-1 + dc.n+l 



), (7) 



which can be identified with the coupled mode equations 
describing optical Bloch oscillations observed in a waveg- 
uide array with linearly growing effective index 
Thus, the recurrences in the GCQW dynamics can be 
understood as a discrete version of o ptic al Bloch oscilla- 
tions. The exact solution of Q (see [171) gives 



m = 4 



[df . 




sin 




- 7T 


\ f J. 



(8) 



Thus, a prefect recurrence should be observed for t = kT. 
In a discrete version, for t = kf and t <C t(p), we indeed 
observe a (nearly) perfect recurrence provided kT is an 
even integer (for each odd t we have Pit) = 0); if this is 
not the case, then an imperfect recurrence occurs at one 



of the neighboring steps. A numerical simulation fully 
confirms these anticipations. 

Finally, let us comment on the case when T = 2ir/§ is 
irrational. Then, T = lim^a, with p n (and q n ) tending 
to infinity as n — > oo. As we have shown above the 
ballistic diffusion can take place only after time r(p n ), 
which in this case goes to infinity, so the diffusion will 
be suppressed forever. Moreover, we can expect that the 
probability that the particle returns to the initial state in 
time t oscillates irregularly with period T. One can also 
easily see that the maximum value cr max of the standard 
deviation (localization length) can be approximated by 



- vl -D. 



(9) 



Again, this agrees well with earlier simulations for D = 
1/2 and f = 4tt given by Buerschaper and Burnett [l^ 
(see Fig. 2d therein, for which a max = 3.4). 

In the paper we have studied the behavior of GCQW, 
where the phase is modified by $ = 2ir^ at each mode. 
For such a system the probabilities of recurrences and 
multiple recurrences are given by PJ and (Jlj respectively. 
They depend strongly on the parity of p which can be at- 
tributed to a "degeneration" of quasi-energy levels in the 
case when p is even. We showed that the system oscil- 
lates for small times t while for t much larger than r(p), 
(see (0), a ballistic diffusion occurs for every D < 1, 
with the standard deviation a given by 10 • For an irra- 
tional <I>/27r, our model predicts dynamical localization 
and gives a good upper bound for the localization length 
Cmax (see ©)• Our results are in perfect agreement with 
experimental and numerical observations from |6lll2| but 
some of their consequences are still to be verified exper- 
imentally (a natural candidate for an experimental real- 
ization of GCQW would be a version of spin-dependent 
transport of atoms in optical lattices |19l l20j ) . 
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